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FUN3D Navier-Stokes solutions were computed for the 4th AIAA Drag Prediction Workshop grid con- 
vergence study, downwash study, and Reynolds number study on a set of node-based mixed-element grids. 
All of the baseline tetrahedral grids were generated with the VGRID (developmental) advancing-layer and 
advancing-front grid generation software package following the gridding guidelines developed for the work- 
shop. With maximum grid sizes exceeding 100 million nodes, the grid convergence study was particularly 
challenging for the node-based unstructured grid generators and flow solvers. At the time of the workshop, 
the super-fine grid with 105 million nodes and 600 million elements was the largest grid known to have been 
generated using VGRID. FUN3D Version 11.0 has a completely new pre- and post-processing paradigm that 
has been incorporated directly into the solver and functions entirely in a parallel, distributed memory envi- 
ronment. This feature allowed for practical pre-processing and solution times on the largest unstructured-grid 
size requested for the workshop. For the constant-lift grid convergence case, the convergence of total drag is 
approximately second-order on the finest three grids. The variation in total drag between the finest two grids 
is only 2 counts. At the finest grid levels, only small variations in wing and tail pressure distributions are seen 
with grid refinement. Similarly, a small wing side-of-body separation also shows little variation at the finest 
grid levels. Overall, the FUN3D results compare well with the structured-grid code CFL3D. The FUN3D down- 
wash study and Reynolds number study results compare well with the range of results shown in the workshop 
presentations. 


I. Introduction 

The fourth international AIAA Drag Prediction Workshop (DPW-IV) was held in San Antonio, Texas in June 2009 
to evaluate CFD transonic cruise drag predictions for subsonic transport aircraft. 1 The objectives of the workshop were 
as follows: 1) to build on the success of past AIAA Drag Prediction Workshops, 2 ^ 1 2) to assess the state-of-the-art 
computational methods as practical aerodynamic tools for aircraft force and moment prediction of industry relevant 
geometries, 3) to provide an impartial forum for evaluating the effectiveness of existing computer codes and modeling 
techniques using Navier-Stokes solvers, and 4) to identify areas needing additional research and development. The 
focus of this particular workshop was on “blind” drag prediction (a priori experimental data was not available for 
comparison) for the NASA Common Research Model (CRM) transonic wing-body-tail configuration. 5 As with prior 
Drag Prediction Workshops, grid convergence of the force and moments coefficients was evaluated at the design cruise 
Mach number and lift coefficient. With maximum grid sizes exceeding 100 million nodes (600 million cells), the grid 
convergence study was particularly challenging for the node-based unstructured grid generators and flow solvers. This 
maximum grid size was four times larger than required for DPW-III. Additionally in this workshop, the drag increment 
due to the trimming of the aircraft and drag increment due to Reynolds number scaling were evaluated. An optional 
Mach sweep study was also requested to look at drag rise. 
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The experimental data for the CRM and the statistical analysis of the workshop results are scheduled for publication 
in June 2010. 6,7 As with previous DPW workshop, the DPW-IV test cases will in all likelihood become standard CFD 
test cases. Eliasson and Peng have already published a paper looking at the influence of turbulence modeling and grid 
resolution in computations of the CRM. 8 Hashimoto et al. have used the CRM to evaluate a hexahedra grid generation 
method. 9 Many more computational studies of the DPW-IV test cases are scheduled for publication at the 28th AIAA 
Applied Aerodynamics Conference in June 2010. 

Results from the DPW-IV grid convergence study, downwash study, and Reynolds number study were submitted to 
the workshop for the FUNS D 10 12 unstructured-grid Navier-Stokes flow solver on a set of node -based mixed-element 
grids. This paper gives a summary of the FUN3D results submitted to the workshop. A description of the unstructured 
node-based grids is provided in this paper detailing the grid parameters and spacings. Some improvements to the 
FUN3D solver/grid processing efficiency and parallelization are also discussed. These improvements allowed for the 
practical processing of the largest unstructured-grid sizes requested for the workshop (100 million nodes). Although 
experimental data are not yet available, comparisons are made with overset structured grid results from the well-known 
code CFF3D. 13 


II. DPW-IV Workshop Cases 

All of the DPW-IV workshop cases are based on the analysis of the NASA CRM wing/body/tail configuration 
shown in Figure 1 . The configuration is a transonic transport aircraft designed to fly at a cruise Mach number of = 
0.85 with a nominal lift condition of Cl = 0.500 at a Reynolds number of Re c = 40 x 10 6 per mean aerodynamic 
chord (c re f). The reference geometric quantities for the CRM are summarized in Table II. A detailed description of 
the CRM can be found in Ref. 5. The first test case (Case 1) for the workshop included a grid convergence study and 
a downwash study. For the grid convergence study, computations were performed for a constant design-lift coefficient 
of Cl = 0.500 ± 0.001 at a freestream Mach number of = 0.85 and freestream chord Reynolds number of 
Re c = 5 x 10 6 . (Note that this Reynolds number represented a planned wind tunnel test condition.) For the grid 
convergence study, a nominal tail incidence angle i/, = 0° was specified. A family of four grids was required for the 
study with consistent levels of refinement between each grid: coarse (3.5 million nodes), medium (10 million nodes), 
fine (35 million nodes), and super-fine (100 million nodes). The downwash study required complete angle-of-attack 
sweeps for the CRM with three different tail settings (/'/, = — 2° , 0° , +2°) on a set of medium sized grids. The trimmed 
drag polar was then derived from an interpolation of the polars at the variable tail incidences using a process provided 
by the organizing committee. The delta drag polar was then computed by differencing the polar results from a tail-off 
configuration with results from the trimmed tail-on drag polar. The optional Reynolds number study (Case 3) is a 
comparison of Re c = 5 x 10 6 results with Re c - 20 x 10 6 results on the medium grid size at the design Mach number 
and lift condition (M^ = 0.85 and Cl = 0.500 ± 0.001). For Case 3, the nominal tail incidence angle ih = 0° was 
used. FUN3D results were not computed for the optional Mach sweep study (Case 2). 


Table 1. Reference geometry for the CRM. 


Mean Aerodynamic Chord ( c re f ) 

275.80 inch 

Wing Reference Area/2 

297, 360 inch 2 

Wing Span/2 

1159.0 inch 

Aspect Ratio! ,4 //) 

9.0 


III. DPW-IV Unstructured LaRC Node-Based Grids 

The FUN3D solutions were computed using a mixed-element version of the vertex-based unstructured tetrahedral 
grids generated for the workshop by Dr. Pirzadeh. The tetrahedral grids are available on the DPW-IV web site under 
the label “unstructured LaRC/NodeBase.” 1 All of the baseline tetrahedral grids were generated with the VGRID (devel- 
opmental) advancing-layer and advancing-front grid generation software package 14 following the gridding guidelines 
developed for the workshop. 1 The grids generated with VGRID were fully tetrahedral. However, VGRID uses an 
advancing layer technique to generate the boundary layer portion of the grid so that prisms can be reconstructed in 
the boundary layer for use with a mixed-element discretization. In the boundary layer, three tetrahedral cells are com- 
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bined into one prism. The mixed-element grids have essentially the same number of nodes and nodal spacing as the 
fully tetrahedral grids although the number of cells and the shape of the control volumes differ in the boundary layers. 
The merging process adds a trivial number of node points in the interface region between the advancing layers and 
advancing front to construct a valid topology in the interface region between the prismatic layers and the tetrahedral 
region. For the LaRC node-based grids, the merging process also reduces the total number of cells in the grid by a 
factor of approximately 1/2 to 2/3. 

For the grid convergence study, a set of four globally-refined grids were generated following the gridding guide- 
lines. Similarly, for the downwash study, a set of three consistently sized grids with varying tail incidences were 
generated. For the tail-off configuration, the tail solid surfaces were removed, and the volume grid was re-generated 
with the same spacing requirements in the location of the tail as were used in the tail-on configuration. A summary 
of the Case 1 mixed-element grid sizes and runs is shown in Table III. For the Case 3 high Reynolds number grid, 
the surface grid distribution from the Case 1 medium grid was used with a smaller minimum wall spacing. The total 
number of grid nodes for the Re c = 20 x 10 6 grid was 1 1 million. Note that at the time of the workshop, the Case 
1 super-fine grid with 105 million nodes and 600 million elements was the largest grid known to have been generated 
using VGRID. 


Table 2. Summary of DPW-IV Case 1 mesh sizes and runs at Moo = 0.85, Re, = 5 x 10®. 


Grid 

Tail Incidence (deg.) 

Grid Size (million nodes) 

C L = 0.500 

Alpha Polar 

Coarse 

0 

3.7 

required 

not required 

Medium 

0 

10.3 

required 

required 

Fine 

0 

35.5 

required 

not required 

Super Fine 

0 

105 

required 

not required 

Medium 

-2 

10.3 

not required 

required 

Medium 

+2 

10.3 

not required 

required 

Medium 

No Tail 

8.2 

not required 

required 


VGRID has two types of spacing requirements: the “inviscid” spacing distributions are used in the advancing- 
front region of the mesh, and the “viscous” spacing distributions are used in the advancing-layer regions of the mesh 
where high stretching is required. The different grids for the grid refinement study were generated by a global coars- 
ening/refinement of the inviscid spacing parameters (VGRID “sources”) and a global coarsening/refinement of the 
viscous wall spacing. A summary of characteristic grid spacing for the wing/body/tail ( ih = 0°) configuration is 
shown in Table 3 to illustrate that the grids were generated following the gridding guidelines for the workshop. The 
global spacing factor shown in Table 3 represents the relative coarsening/refinement of the inviscid spacing parameters 
for each of the grids. The range of chordwise spacing for the wing and tail leading edge (LE) and trailing edge (TE) 
fall close to the 0.1% local chord specified for the medium grids at the root chord ( c r0 ot ) and tip chord ( cu p ) locations. 
The spanwise stretching of the grids was 10:1 near the wing LE and 15:1 near the wing TE. There were similar limits 
on the spanwise stretching near the tail LE and TE. The limits on the grid spanwise stretching rates keep the wing 
and tail spanwise spacing at root and tip well within the workshop requirement (0.1% local semi-span for the medium 
grid). A composite view of the medium grid wing/body/tail grid is shown in Figure 2 to illustrate the overall topology 
of the surface mesh. 

The clustering of points normal to the surface was computed according to the VGRID stretching function 15 

S n = S l [l + r 1 (l + r 2 ) n - 1 ] n - 1 (1) 

where S n is the normal spacing of the n th layer, <5i is the spacing of the first layer, and the factors r ± and r -2 are constants 
that determine the rate of stretching. (Note if r 2 is zero the stretching is geometric.) Based on the advancing layer 
stretching factors given in Table 3, the maximum growth rate of the cell sizes in the viscous layers is approximately 
1.19 which is less than the 1.25 requested in the gridding guidelines. The blunt trailing edges of the wing and nacelles 
were resolved in all grids by explicitly dividing the trailing-edge surfaces into the specified number of cells. 
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Table 3. Summary of characteristic grid spacing for the DPW-IV Unstructured LaRC Node-Based wing/body/tail if, = 0° grids. 




Node-Centered Grids 



Coarse 

Medium 

Fine 

Super Fine 

Global Spacing Factor 

1.654 

1.00 

0.562 

0.355 

% Chordwise Spacing 

0.20 % Croat 

0.12% c root 

0.01%c root 

0.04% c root 

at wing LE 

0.20 %ct ip 

0.l2%c Up 

0.07 %c t i P 

0. 04% c t jp 

% Chordwise Spacing 

0.07% Croat 

0.04 %c root 

0.02 %c root 

0.01% c root 

at wing TE 

0.12 %c t ip 

0.01%c Up 

0.04 %c tip 

0.02 %c tip 

% Chordwise Spacing 

0.20% Croat 

0. 16 %c roo t 

0.09 %c root 

0.06 %C r0 ot 

at tail LE 

0. 26% Cup 

0.12 %c tip 

0.07 %c t ip 

0.04 %c tip 

% Chordwise Spacing 

0. 20% C root 

0.12 %c root 

0.07%c root 

0.04 % Croot 

at tail TE 

0. 20% c t i P 

0.l2%c t i P 

0.07 %c tip 

0.04 %ctip 

Wing/Tail TE Cells 

8/8 

12/12 

16/16 

24/24 

Spacing Fuselage Nose/Aft Body 

4.0% C re y 

2.4% c re / 

1.3% c re / 

0.9 %c ref 

Average Cell y + 

1.0 

0.667 

0.444 

0.296 

Nominal Prism 
Layer Cells 

25 

25 

25 

25 

Viscous Wall 
Spacing <5y (inches) 

0.001478 

0.0009850 

0.0006570 

0.0004380 

Viscous Stretching 

0.12, 

0.12, 

0.12, 

0.12, 

rr,r 2 

0.02 

0.02 

0.02 

0.02 

Outer Boundary Box 

100 c re f 

100 c re f 

100 c re / 

100 c re f 


IV. FUN3D Flow Solver 

FUN3D 111 12 is a finite-volume Reynolds-averaged Navier-Stokes (RANS) solver in which the flow variables are 
stored at the vertices or nodes of the mesh. FUN3D solves the equations on mixed element grids, including tetrahedra, 
pyramids, prisms, and hexahedra and also has a two-dimensional path for triangular/quadrilateral grids. It employs 
an implicit upwind algorithm in which the inviscid fluxes are obtained with a flux-difference-splitting scheme. At 
interfaces delimiting neighboring control volumes, the inviscid fluxes are computed using an approximate Riemann 
solver based on the values on either side of the interface. Several convective flux schemes are available in FUN3D. 
The most common scheme for subsonic and transonic flows is Roe’s flux difference splitting 16 which is used in 
the current study. For second-order accuracy, interface values are obtained by extrapolation of the control volume 
centroidal values, based on gradients computed at the mesh vertices using an unweighted least-squares technique. 
Several flux limiters are available in FUN3D to limit the reconstructed values when necessary. The most common used 
for transonic flows is that of Venkatakrishnan 17 which is used in this study. For tetrahedral meshes, the full viscous 
fluxes are discretized using a finite-volume formulation in which the required velocity gradients on the dual faces are 
computed using the Green-Gauss theorem. On tetrahedral meshes this is equivalent to a Galerkin type approximation. 
For non-tetrahedral meshes, the same Green-Gauss approach can lead to odd-even decoupling. A pure edge-based 
approach can be used to circumvent the odd-even decoupling issue, but yields only approximate viscous terms. Thus 
for non-tetrahedral meshes, the edge-based gradients are combined with Green-Gauss gradients, which improves the 
h-ellipticity of the operator, and allows the complete viscous stresses to be evaluated. 18, 19 This formulation results in 
a discretization of the full Navier-Stokes equations. 

The solution at each time-step is updated with a backwards Euler time-differencing scheme. At each time step, the 
linear system of equations is approximately solved with either a multi-color point-implicit procedure or an implicit- 
line relaxation scheme. 20 Local time-step scaling is employed to accelerate convergence to steady-state. For turbu- 
lent flows, several models are available within FUN3D. The most common are the one-equation model of Spalart- 
Allmaras 21 (S-A) and the two-equation shear stress transport (SST) model of Menter. 22 The S-A model may be solved 
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loosely coupled to the mean-flow equations or tightly coupled to the mean-flow equations. The Menter SST Model 22 
is loosely coupled. The loosely coupled S-A model is used for the current study. 

Efforts were initiated in the summer of 2006 to improve the computational performance of the FUN3D code. 
Many low-level aspects of the code were examined including cache reuse, MPI communication, inlining, basic blocks 
and alternative ordering techniques for the grid/linear algebra operations. The result of this effort was a factor of 6.5 
speed-up of the recently-released Version 11.0 over Version 10.3.1 (May 2006). The code performance at several 
points in time from May 2006 to the the current Version 11.0 is shown in Figure 3(a). The comparison of density 
residual versus wall clock time shows a 6.5 times speed-up for a wing/body transonic Navier-Stokes solution to the 
same level of residual convergence. This timing study was performed with fixed compiler, compiler options, computer 
hardware and solver options so that all of the code speed-up shown is due to computer science and implementation 
differences. The FUN3D solver has also demonstrated linear scaling in parallel efficiency up to 8,192 computer cores 
on the Pleiades supercomputer at the NASA Advanced Supercomputing (NAS) facility. Pleiades is an SGI ICE system 
with Xeon E5472 and X5570 processor cores. The scaling study utilized the Xeon E5472 processors with 2 quad-core 
processors per node and 8 GB of memory per node. Figure 3(b) shows the parallel scalability of the Version 1 1.0 code 
on a range of grid sizes including the DPW-IV 105 million node super- fine mesh described in the previous section. 
Grids from the second DPW (DPW-II) and the Ares- IX launch vehicle are also included in Fig. 3(b). Most of the 
results shown were run fully-dense; namely, one process per core (e.g., eight processes on a dual-quad node). One 
timing study for the Ares- IX was run in a half-dense environment (four processes on a dual -quad node) to illustrate 
the speed-up from decreasing the inter-processor communication. The solutions for the the timing study were not run 
to convergence so no absolute times are provided in Fig. 3(b). 

FUN3D Version 11.0 also has a completely new pre- and post-processing paradigm that has been incorporated 
directly into the solver and functions entirely in a parallel, distributed memory environment. The result is an end-to- 
end scalable approach that eliminates hardware constraints as problem sizes grow. The new paradigm also cuts down 
on pre- and post-processing times considerably and greatly simplifies the use of FUN3D. For the 105 million node 
tetrahedral grid, grid partitioning and pre-processing took approximately 8 minutes using 1,024 cores on Pleiades. 
Partitioning the same grid on 1,024 processors using serial pre-processor previously required two weeks using 800 
gigabytes (GB) of shared-memory on the NAS Columbia supercomputer which has SGI Altix 3700/4700 cores with 
approximately 2GB of memory per core. A substantial decrease in I/O processing time, and I/O memory overhead, 
was seen using either of the two new I/O options within Version 1 1.0 (i.e., FORTRAN stream I/O, MPI/IO). 

V. Computational Results 

The FUN3D workshop results were computed using the full Navier-Stokes equations with the flux-difference- 
splitting scheme of Roe, the Venkatakrishnan flux-limiter, and the loosely-coupled S-A turbulence model on mixed 
element grids (prisms through the boundary layer, tetrahedra in the farfield). In accordance with the workshop spec- 
ifications, all computations assume fully turbulent flow. For all cases, local time-step scaling was used to accelerate 
the residual convergence of the non-linear equations and achieve steady-state values of forces and moments. For the 
constant lift cases, the angle of attack {AO A) was relaxed based on the difference in total lift and the target lift during 
the convergence history. (Note that updating the angle of attack at every iteration slows the solution convergence 
compared to a case where the angle of attack is constant.) Good solution residual, force, and moment convergence 
was achieved for the all the workshop computations even on the super-fine grid size. Figure 4 shows the non-linear 
solution residual, force, and moment convergence histories for the medium and super-fine grids. The F2 norm of the 
density (R±) residual and the turbulence model (Rq) residual are shown for each local time-step. 

Total force and moment results for the CRM grid convergence study are given in Fig. 5. The angle of attack, total 
drag ( Cd X pressure drag ( Cd p ), skin friction drag ( Cd s /) and pitching moment (Cm?/) are plotted versus N ~ 2 / 3 , 
where N is the number of nodes for the unstructured grid. (In the asymptotic range, one would expect a linear 
variation in forces or moments with N~ 2 / 3 for a second order scheme for a family of three-dimensional grids that 
are globally refined.) Thus, results using finer grids appear to the left in the figures, and results using coarser grids 
appear to the right. The variation in total drag on the finest three grids is about 6 counts, and the variation on the 
finest two grids is 2 counts. The decomposition of the drag into pressure and skin friction components illustrates that 
the majority of the variation is due to the pressure component. While the changes in drag and angle of attack are 
nearly linear on the finest three grids, the change in pitching moment is not, possibly due to the sensitivity of the wing 
pitching moment to shock strength. Since the experimental data and workshop statistics are not currently available, 
the FUN3D results from the grid convergence study are compared to results from a well-known structured-grid code, 
CFF3D. 23 The CFF3D results were computed on a set of overset grids using an upwind algorithm, the SA turbulence 
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model, and the thin-layer approximation of the RANS equations. A complete description of the overset grids and 
CFL3D results can be found in Ref. 13 . The results for the CFL3D DPW-IV grid convergence study are compared 
with the FUN3D results in Fig. 5. Overall, the comparison between the computed forces and moments is good. The 
comparison in total drag between the CFL3D results on finest grid (55M cells) and the FUN3D results on the fine and 
super-fine grids is within 5 drag counts with most of the variation coming from the pressure component. Similarly, the 
comparison in pitching moment is quite good at the finest grid levels. 

The component forces and moments for the wing, body and tail are shown in Fig. 6. The CFL3D results from 
Ref. 13 are also included for comparison. The range of the vertical scales is the same for each component of lift, drag 
and pitching moment which allows for a comparison of variation between the component force and moment contribu- 
tions. (The data ranges for drag and pitching moment in Fig. 6 are also the same as shown in Fig. 5.) Figures 6(a)-(c) 
show that the lift contributions from each of the components is relatively constant with grid refinement. The drag on 
each of the components consistently decreases with grid refinement as shown in Figures 6(d)-(f). However, the tail 
variation is much less than the wing and body variation which are similar in magnitude. Figures 6(g)— (i) show that 
the pitching moment contributions from the body and tail are relatively constant with grid refinement with most of the 
variation coming from the wing. In comparing the FUN3D results to the CFL3D results, the structured grid solutions 
consistently have more of a download on the tail which explains the prediction of a more nose-down tail pitching 
moment from CFL3D. Figures 6(d)— (f) illustrate that the body and tail drag predictions are very close on the fine grids 
with most of the difference in total drag between the FUN3D and CFL3D results come from the wing component. 

Figures 7 and 8 show the grid convergence of the wing and tail chordwise pressure coefficient ( C p ) distributions, 
respectively. The computational results are shown at selected stations across the wing and tail span that correspond 
to requested workshop and experimental data locations. (Note that the inset upper-surface pressure line contours are 
from the fine grid solution.) Figure 7 shows that on the inboard half of the wing the differences between the surface 
pressures on all the grid densities is very small. On the outboard wing, there are more differences in shock strength 
across the grid range, especially at the tip where the “double shock” structure does not show up at the coarsest grid 
density. However, the pressure distributions on the fine and super fine grids are very close. Figure 8 shows that 
the variation in surface pressure with grid refinement is very small across the tail with only small differences in the 
leading edge pressure peak. The wing and tail surface pressures from the FUN3D medium and fine grid solutions 
are compared to the CFL3D medium grid solutions in Fig. 9. (Note that the fine and medium grid CFL3D pressures 
distributions are very similar to within plotting accuracy, so only the medium grid results are shown.) The wing 
pressures at the mid-span location compare very well between the two codes. At the wing tip location, the overall 
pressure distribution compares well except near the double-shock feature, where FUN3D’s two shocks are positioned 
slightly closer together. The tail pressures at the mid-span and tip locations also compare very well between the two 
codes with only a slight decrease in nose pressure on the lower surface of the tail for the CFL3D results. This difference 
could explain the increased tail load predicted with CFL3D. 

The CRM fine grid skin friction coefficient and surface restricted streamlines shown in Figure 10 illustrate the 
flow patterns on the wing and tail upper surface. The flow patterns on the fine grid are indicative of those seen on the 
other grid densities. As was seen in DPW-II and DPW-III, there is an area of flow separation at the trailing edge of the 
wing-root juncture. There is also a pattern of small trailing-edge separation on the outboard portion of the wing. The 
grid convergence of CRM wing chordwise skin friction distribution at the r] = 0.603 span station shown in Figure 1 1 
illustrates that there is a small increase in trailing-edge separation with grid refinement. The effect of grid refinement 
on the trailing-edge wing-root juncture flow is shown in Figure 12. In this figure, a close-up of the surface restricted 
streamlines on the wing upper surface is shown from a perspective just downstream and outboard of the trailing edge. 
The comparison of streamline patterns shows an increase in the size of the separated region with grid refinement 
between coarse, medium, and fine grids, but then very little change between fine and super-fine. No trailing-edge or 
side-of-body separation was noted on the tail at any grid density. Not all codes in the workshop predicted the same 
flow patterns on the CRM. 1 Although some codes in the workshop predicted similar patterns of flow separation on 
the CRM, other codes predicted no wing-root juncture separation and/or no wing trailing-edge separation. Also some 
codes predicted a separation pattern on the fuselage just downstream of the tail trailing edge. Although not shown, 
CFL3D predicted a slightly smaller side-of-body separation region than FUN3D at the medium grid resolution, but 
the bubble decreased in size as the grid was refined. This may be due to the use of the thin-layer approximation in the 
CFL3D results. Use of the thin-layer approximation has been shown to reduce the size of the side-of-body separation 
in DPW-III computations. 24 

Figures 13 and 14 show the computational results of the downwash study. Figure 13 illustrates the effects of the 
tail incidence on the CRM total lift, idealized drag, and pitching moment coefficients, where idealized drag is defined 
as Ci;> — C\j (~ AH). The CFL3D results from Ref. 13 are also included in Fig. 13 for comparison. The comparison of 
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lift prediction between the two codes is quite good across the angle-of-attack and tail incidence range. Similarly, the 
comparison of pitching moment is also quite good except at the highest angle of attack. The idealized drag prediction 
from CFL3D is consistently lower than the FUN3D results across the angle-of-attack and tail incidence range except 
at 4 degrees angle of attack where they are very close. It is interesting to note that even when the codes do not predict 
the same absolute values of forces and moments the deltas due to tail incidence appear to be at the same levels. The 
trimmed polar forces and moments interpolated from the FUN3D results are shown in Figure 14 along with the tail-off 
polar data. The delta drag due to trimming the tail at the design lift coefficient Cl = 0.500 is 27 counts. Range of 
delta drag due to trimming from the workshop presentations was 27-40 counts. 1 

For the Case 3 Reynolds number study as the Reynolds number was increased to Re c = 20 x 10 6 , the angle of 
attack decreased from 2.32° to 2.14°, and the total CRM drag dropped by 33 counts. The drag due to pressure dropped 
by 21 counts, and the drag due to skin friction dropped by 12 counts. The wing-root juncture separation decreased by 
about half the size. The range of delta drag due to increased Reynolds number from the workshop presentations was 
30-35 counts. 1 


VI. Summary 

FUN3D Navier-Stokes solutions were successfully computed for the DPW-IV grid convergence study, downwash 
study, and Reynolds number study on a set of node-based mixed-element grids. All of the baseline tetrahedral grids 
were generated with the VGRID (developmental) advancing-layer and advancing-front grid generation software pack- 
age. Although the grids generated with VGRID were fully tetrahedral, a pre-processing step was used to combine 
the tetrahedral cells in the boundary-layer portion of the grid into prisms. For the grid convergence study, a set of 
four globally refined grids were generated following the gridding guidelines. Similarly, for the downwash study, a 
set of three consistently sized grids with varying tail incidences were generated. For the tail-off configuration, the 
tail solid surfaces were removed, and the volume grid was re-generated with the same spacing requirements in the 
location of the tail as were used in the tail-on configuration. For the high Reynolds number grid, the surface grid 
distribution from the downwash study medium grid was used with a smaller minimum wall spacing. With maximum 
grid sizes exceeding 100 million nodes, the grid convergence study was particularly challenging for the node-based 
unstructured grid generators and flow solvers. At the time of the workshop, the super-fine grid with 105 million nodes 
and 600 million elements was the largest grid known to have been generated using VGRID. Some improvements to the 
FUN3D solver/grid processing efficiency and parallelization were shown which allowed for practical pre-processing 
and solution times on the largest unstructured-grid size requested for the workshop. 

The FUN3D workshop results were computed using the full Navier-Stokes equations with the flux-difference- 
splitting scheme of Roe, the Venkatakrishnan flux-limiter, and the loosely-coupled S-A turbulence model. Good 
residual, force, and moment convergence was achieved for the all the workshop computations even on the super-fine 
grid size. Overall, the FUN3D results compared well with the structured-grid code CFL3D on the finer grid levels. 
For the constant lift grid convergence case, the convergence of total drag was approximately second-order on the finest 
three grids. The variation in total drag on the finest three grids was about 6 counts, and the variation on the finest 
two grids was 2 counts. At the finest grid levels, only small variations in wing and tail pressure distributions were 
seen with grid refinement. Similarly, a small wing/side-of-body separation also showed little variation at the finest 
grid levels. The FUN3D downwash study and Reynolds number study results compared well with the range of results 
shown in the DPW-IV workshop presentations. Additional comparisons of the FUN3D results will be made when the 
experimental data for the CRM and the statistical analysis of the workshop results are made publicly available. 
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Figure 1. NASA Common Research Model wing/body/tail configuration. 
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Figure 2. Surface grid for the CRM wing/body/tail medium-sized grid (f 0.3 million nodes). 
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(a) Computational Timing 



Number of Cores 


(b) Parallel scalability 

Figure 3. FUN3D computational performance. 
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(a) Medium grid (10.3 million nodes) 


Const. AOA 




(b) Super-Fine grid (105 million nodes) 

Figure 4. Non-linear solution residual, force, and moment convergence history for the CRM wing/body/tail. 
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Figure 5. Grid convergence of CRM total forces and moment at Moo = 0.85, Re c = 5 x 10 6 . 
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(g) Wing Pitching Moment (h) Body Pitching Moment (i) Tail Pitching Moment 


Figure 6. Grid convergence of CRM component forces and moment at Moo = 0.85, Re c = 5 x 10 6 . 
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Figure 7. Grid convergence of CRM wing chordwise pressure distributions at = 0.85, Re, = 5 x 10 6 (surface pressure contours 
from fine grid solution). 


14 of 21 


American Institute of Aeronautics and Astronautics 







x/c x/c 



Figure 8. Grid convergence of CRM tail chordwise pressure distributions at Moo = 0.85, Re c = 5 x 10 6 (surface pressure contours from 
tine grid solution). 


15 of 21 


American Institute of Aeronautics and Astronautics 


Medium Grid (10M) FUN3D 
Fine Grid (35.0M) FUN3D 
Medium Overset Grid (16M) CFL3D 



0.00 0.20 0.40 0.60 0.80 1.00 



x/c 


x/c 



Figure 9. Comparison of FUN3D and CFL3D CRM wing/tail chordwise pressure distributions at Moo = 0.85, Re c = 5 x 10®. 
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Figure 10. CRM fine grid (35 million nodes) skin friction coefficient and surface restricted streamlines at M a 0 = 0.85, Re c = 5 x 10 6 . 
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Figure 11. Grid convergence of CRM wing chordwise skin friction distribution at M oo = 0.85, Re, = 5 x 10®. 
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(c) Fine Grid 



(b) Medium Grid 



(d) Extra-Fine Grid 


Figure 12. Effect of grid density on the CRM {Cl = 0.5) trailing edge wing-root juncture separation. 
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Figure 13. Effect of tail incidence on CRM forces and moments at M = 0.85, Re c = 5 x 10 6 . 
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Figure 14. Comparison of the CRM tail-off and trimmed tail-on forces and moments at Moo = 0.85, Re c = 5 x 10®. 
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